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Recent developments in lattice QCD calculation of flavor singlet nucleon matrix elements are reviewed. Sub- 
stantial sea quark contributions are found in the tt-N a term and the quark spin content of the nucleon such that 
the total magnitude including valence contributions is in reasonable agreement with experiments. Some problems 
with flavor non-singlet nucleon matrix elements are pointed out. Recent work on lattice QCD calculation of 
hadron scattering length is also discussed. 



1. INTRODUCTION 

While extensive effort has been invested for 
studies of the flavor non-singlet sector of strong 
interactions with the use of lattice QCD, attempts 
toward understanding the flavor singlet sector 
were quite limited until a few years ago . This 
stemmed from the severe technical difficulty that 
a calculation of disconnected quark loop ampli- 
tudes, not present in the flavor non-singlet sec- 
tor, requires quark propagators for a number of 
source points equal to the space-time lattice vol- 
ume, which is prohibitively computer time con- 
suming if it is to be carried out with the conven- 
tional method of point source. 

Several methods have been developed in the 
last few years, however, for overcoming this prob- 
lem. One of the method is a variant version of the 
wall source technique in which gauge configura- 
tions are not fixed to any gauge Q, as was em- 
ployed in the early studies of extended sources . 
Another is an improvement of the method of noisy 
source with the use of a random Z2 noise 

The first method has been successfully applied 
to estimate the flavor singlet r]' meson mass in 
quenched QCD The method also allows a 
calculation of full hadron four-point functions, for 
which a similar computational problem is present, 
and hence an evaluation of hadron scattering 
lengths These results were available at 

the time of the Dallas Conference |ic| ]. 

In this article we review new developments in 
this field since then. We devote the major part of 



the review to discuss recent calculations of flavor 
singlet nucleon matrix elements [|ri]-|l^,?], specif- 
ically the TT-N a term and the axial vector matrix 
elements (quark content of proton spin) . We also 
point out some problems with flavor non-singlet 
nucleon matrix elements which were not ap- 
parent in the pioneering studies pO| , pT[ |. In the 
rest of the article we turn to some recent work on 
hadron scattering lengths . 

Most of calculations discussed in this review 
are carried out in quenched QCD with the Wilson 
quark action. This should be understood unless 
stated otherwise. 

2. CALCULATIONAL TECHNIQUES 

A standard method for extracting the nucleon 
matrix element of an operator O is to employ the 
formula for the ratio of three- to two-point func- 
tion of the nucleon and the operator given by |20| 



(^WE„oN^(o)) 

{N{t)Nm 

const + Z^^{N\0\N) t, 



(1) 



where N(t) is the nucleon operator projected to 
zero momentum and Zq is the lattice renormal- 
ization factor for the operator O. For a flavor 
singlet quark bilinear operator O = qTq, there 
are two types of diagrams contributing to the 
three-point function shown in Fig. |^. The con- 
nected amplitude can be calculated by the well 
known source technique psl. To handle the dis- 
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Figure 1. (a) Connected and (b) disconnected 
contribution to the three-point function of nu- 
cleon and flavor singlet operator qVq. 

connected piece, for which the source point of the 
disconnected quark loop has to be summed over 
all space-time sites, two methods have been de- 
veloped. 

2.1. Wall source method without gauge fix- 
ing 

In this method one calculates quark propaga- 
tors with unit source at every space-time site 
without gauge fixing p|,pl|. The result G{n) is 
a sum of point-to-point propagators of form 



G(n) = ^G(r 



(2) 



The product of the nucleon propagator and 
J2n Tr[G(n)r] equals the disconnected amphtude 
up to gauge-variant nonlocal terms, which, how- 
ever, cancel in the average over gauge configu- 
rations. Generalizing the well-known argument 
, it is possible to show that the magnitude of 
gauge- variant noise terms is smaller than the sig- 
nal by a factor 1/^/VN for a sufficiently large 
volume V and a number of configurations TV. 
In practice 0(100 — 200) configurations are em- 
ployed to obtain reliable signals on a lattice such 
as 16^ X 20 at /? = 5.7 ]l2|,^. 

This method requires only two quark matrix in- 
versions for each gauge configuration, one for cal- 
culating the nucleon propagator and the other for 
the disconnected quark loop. Repeating the latter 
n times, applying random gauge transformations 
to a given gauge field, does not yield much gain 
since errors are reduced only by a factor l/y/n 
while the computer time grows proportional to n. 



2.2. Z2 source method 

In this method ||^ L quark propagators are 
calculated for an L set of random Z2 source 77^ 
placed at every space-time site for each gauge con- 
figuration, i.e., 



(3) 



The product of the nucleon propagator and 
T Tr[Gf(?7^r] equals the disconnected ampli- 
tude since hmi^oo x Y.i=i vivi' = ^nn'- 

Employing a random source for evaluating dis- 
connected amplitudes was previously tried for the 
TT-N a term for the Kogut-Susskind quark ac- 
tion without much success . The improvement 
with the Z2 noise is that the variance of the result 
for a finite L is smallest compared to other types 
of noise including the Gaussian noise. 

The Z2 source method allows a calculation of 
disconnected amplitudes for each configuration. 
However, concrete applications of the method 
show that the number of quark propa- 
gators L has to be large, generally exceeding 
100. Thus the method becomes harder to use for 
lighter quarks for which a substantial computer 
time is necessary to obtain even a single quark 
propagator. 

2.3. Other methods 

Some fiavor singlet nucleon matrix elements 
can be estimated indirectly employing their re- 
lation with other quantities: (i) the scalar den- 
sity matrix element can be related to the nucleon 
mass via Hellmann-Feynman theorem. 



{N\qq\N) = 



drriN 
dm„ 



(4) 



(ii) for the axial vector matrix element the U{1) 
anomaly equation yields 



{P,s\qj^j5q\p, 



(X hm ^— ,(p + q,s\ TtF^^F^^ |p, s) 
q^a q ■ s 



(5) 



Application of these relations require full QCD 
simulations if disconnected contributions are to 
be included. In particular the anomaly relation 
(H) fails in quenched QCD due to the double pole 



Table 1 



singlet matrix elements 

Kyoto-Tsukuba |ll-|ll 
Kentucky M M 



Altmeyer et al. |17| 
Gupta et al. full ||l8| 



/3 


77l7r / nip 


size 




^ conf. 


5.7 (quenched W) 


0.6 - 0.86 


16^ 


X 20 


260 


6.0 (quenched W) 


0.8 - 0.95 


16^ 


X 24 


50 


5.35 {Nf = 4 KS) 


0.46 


16^ 


X 24 


85 


5.4-5.6 {Nf = 2 W) 


0.7 - 0.85 


16^ 


X (16 X 2) 


15 - 45 



non-singlet matrix elements 
Gupta et al. |l|l 6.0 (quenched W) 

Liu et al. [|§ 6.0 (quenched W) 

Gockeler et al. Q 6.0 (quenched W) 



0.72 - 0.79 
0.8 - 0.95 
0.7 - 0.89 



16^ X 40 
16^ X 24 
16^ X 32 



35 

24 

400-1000 



in the rj' meson propagator connecting the discon- 
nected quark loop and the nucleon |19|] . For these 
reasons, only a few attempts have been made with 
these methods in the past [ ^^|jr^ . 

3. NUCLEON SCALAR MATRIX ELE- 
MENTS 

3.1. TT-N <T term phenomenology 

The TT-N a term is defined as the product of 
the nucleon matrix element of the scalar density 
and the average quark mass m = (m„ -f mrf)/2 



a{t) = m{N{'[l)\uu + dd\N{v)) 



(6) 



evaluated at i = {p' —pY = 0, i.e., a = cr(0). Cur- 
rent algebra and PCAC relate a{t) to the crossing 
even ir-N scattering amplitude at the unphysi- 
cal Cheng-Dashen point t = Irn^ and a disper- 
sion analysis of ir-N scattering leads to the value 
a{2ml) = 64(8)MeV [||,|6|. 

On the other hand one can write 



a(0) 



(To m{N\uu + dd'- 2ss\N) (7) 



with 



y 



2{N\ss\N) 
(N\uu + dd\N)' 



(8) 



Treating flavor SU{3) breaking to first order, one 
can estimate [ETl 



CTO 



{Me + Ms - 2Mn) « 25MeV, (9) 



where nis is the strange quark mass. If one as- 
sumes that the variation A^- = a{2'm^) — a{0) is 



small, the above two estimations lead to a large 
strangeness content for the nucleon y w 0.6. 

It has been argued, however, that the variation 
is substantial, A^- « 15 MeV |^^, in which case 
the value of the a term is reduced to cr « 45 MeV. 
Combined with the suggestion |^ that one-loop 
chiral perturbative corrections raise the value of 
(Tq to (To ~ 35MeV, this implies a more reasonable 
value y « 0.2 for the nucleon strangeness content. 

Let us comment that the magnitude of varia- 
tion Ao- may be examined within lattice QCD by 
extrapolating a{t) from the physical region t < 0. 
A first attempt has been reported at this confer- 
ence [|9). 

3.2. TT-N a term recent lattice results 

In Table we list recent studies of nucleon ma- 
trix elements together with parameters of runs. 

Results obtained for the connected contri- 
bution in the scalar matrix element {N\uu -\- 
Jd|iV)conn. are plotted in Fig. ||, where we used 
the tadpole-improved renormalization factor in 
the MS scheme at the scale /i = 1/a given by |32j 



3K 

4i^: 



1. 



l- 0.0098a^(-) 



(10) 



The "experimental" value at {rriT^/mp)'^ = 
equals {N\uu + dd\N) = a/rh calculated with 
cr = 45 MeV and m = 5.5 MeV. Clearly the con- 
nected contribution is too small to account for the 
experimental value of the tt-N a term. 

Previous results for the connected contribution 
pc| , |2l| are consistent with those in Fig. ^. Esti- 
mates from the nucleon mass in quenched QCD 
using the relation (4) also give consistent values. 
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Figure 2. Connected contribution {N\uu + 
dd\N)conn. for the scalar matrix element as a func- 
tion of (rriTt/mpf' . For reference see Table 1. 

A direct calculation of the disconnected ampli- 
tude (Fig. |l|(b)) has been made in quenched QCD 
using both the wall source method without gauge 
fixing ||ll|,|l^ and the Z2 source method Q,^. 
In Fig.j^results for the R ratio (|l|) obtained with 
the wall method ||ll|,|l|l at /? = 5.7 and K = 0.164 
{rriTr/mp = 0.74) are plotted. We observe reason- 
able signals with a linear increase in t for the dis- 
connected amplitude(Fig. §(a)), albeit errors are 
significantly larger compared with those for the 
connected amplitude(Fig. ||(b)). 

The disconnected amplitude obtained with the 
Z2 source at /3 = 6.0 and K = 0.154 

(m-^lmp = 0.80) is shown in Fig. ^. A reason- 
able linear increase is also seen for this case. 

An important outcome of these calculations is 
that the disconnected contribution to the scalar 
matrix element is large, as is clear from a similar 
magnitude of slope in Fig. || (a) and (b). Further- 
more the total magnitude of the matrix element 
obtained by adding the disconnected contribution 
to the connected one is comparable to the phe- 
nomenological estimate. This is shown in Fig. ^ 
where solid symbols represent total values and 
open ones the connected contribution duplicated 
from Fig. H. A linear extrapolation of quenched 
results to the chiral limit yields values of the ma- 
trix element tabulated in Table 2 (a). 

The filled triangles in Fig. ^ represent data 
obtained in two-flavor full QCD using the nu- 
cleon mass derivative dm n /dniq JlSt , showing an 
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Figure 3. Ratio R{t) for the (a) disconnected and 
(b) connected amplitudes of the scalar matrix ele- 
ment {N\uu + dd\N) evaluated by the wall source 
method without gauge fixing for _ft'=0.164 

at (3—5.7 on a 16^ x 20 lattice. Solid lines are 
linear fits over A < t < 9. 
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Figure 4. Ratio R{t) for the disconnected am- 
plitude of the scalar matrix element (N\uu\N) 
evaluated by the Z2 source method for 
K=0.1M at 13^6.0 on a 16^ x 24 lattice. Solid 
line is a linear fit over 8 < i < 12. 



5 



agreement with quenched estimates (fined circles 
and diamonds) in the region of heavy quark with 
m^/mp > 0.7-0.8. 

We find these resuhs to be quite encouraging. 
On a closer inspection, however, we notice several 
problems, to which we now turn. 

3.3. Problems with scalar matrix elements 

3.3.1. physical value of the tt-N a term and 
quark mass 

To convert results for the scalar matrix ele- 
ments into those for the a term, we need the 
value of the quark mass rh in physical units [^3| . 
In quenched QCD current lattice estimates are 
m = 5 — 6MeV for the Wilson quark action, in an 
agreement with the phenomenological estimation 
of rh =5.5MeV. Combined with the total values of 
the scalar matrix element in Table 2 (a), quenched 
lattice results for the a term are consistent with 
the phenomenological value cr « 45 MeV. 

On the other hand, full QCD simulations with 
the Wilson quark action yield rh~ 2 — 3MeV (see 
e.g., Ref. (l^), which are a factor 2 — 3 smaller 
than in the quenched case. Since values of matrix 
elements are similar between the two cases, this 
means that the a term is also small a « 20MeV 
in full QCD. 

Let us add that the a term estimated from the 
nucleon mass through (4) in full QCD with the 
Kogut-Susskind quark action is similarly small 
cr « 20MeV. For this action the quark mass ob- 
tained for quenched and full QCD is consistent, 
being in the range ?fi « 2 — SMeV for both cases. 

Origin of the small value of quark mass and 
that of the a term in full QCD is not clear at 
present. The discrepancy of quark mass between 
quenched and full QCD with the Wilson quark 
action is also not understood. 

3.3.2. strange quark content within nu- 
cleon 

The nucleon matrix element of the strange 
quark density (A^|ss|A^) receives contribution 
only from the disconnected amplitude. In Table 2 
(b) we summarize recent results [^0. These 
values are obtained by extrapolating the valence 
quark mass to the chiral limit while keeping the 
strange quark mass fixed at the physical value. 
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Figure 5. Total amplitude of the scalar matrix el- 
ement {N\uu + dd\N) as a function of (m„/mp)^. 
For reference see Table 1. Triangles arc full QCD 
results obtained with the relation (4). 

Table 2 

Results for nucleon scalar matrix elements 





Ref. y 


Ref. y 


(a) {N\uu + dd\N) 






connected 


2.52(6) 


2.85(6) 


disconnected 


6.2(1.9) 


3.79(13) 


total 


8.9(1.9) 


6.6(2) 


(b) strangeness content 




{N\ss\N) 


2.89(61) 


1.98(9) 


ms{N\ss\N)/mN 


0.40(8) 


0.27(1) 


_ 2{N\ss\N) 
y ~ (N\uu+dd\N) 


0.65(20) 


0.60(4) 



Although the strange quark density within nu- 
cleon is not directly measurable by experiments, 
both the ratio 'ms{N\ss\N)/'mN, which naively 
measures the strange quark contribution to the 
nucleon mass, and the y parameter appear quite 
large. 

Concerning this point Lagae and Liu suggested 
1 34] that the quark mass correction factor 1-f m^a 
1 35], derived for valence quarks nearly on the mass 
shell, may not adequately account for 0{mqa) 
effects for sea quarks in disconnected diagrams, 
which are mostly off the mass-shell states. To 
examine this possibility they calculated the ratio 
Ks{mqa) of the triangle diagram for the opera- 
tor qq with two external gluons at zero momenta 
on the lattice and in the continuum in one-loop 
perturbation theory. 

In Fig. ^ we compare the m^a dependence 



6 



with l+rriga correction 
' with K_ correction 



1.4 
1.2 
1.0 
0.8 
0.6 
0.4 
0.2 
0.0 

0.00 0.05 0.10 0.15 0.20 0.25 0.30 

Figure 6. y parameter as functions of nisa with 
1 + msa or Ks{msa) factor for strange quark den- 
sity. 

of the y parameter for the data of the Kyoto- 
Tsukuba group corrected by the 1 + nisa factor 
and by the Lagae-Liu ks factor. With the ks 
factor the y parameter is reduced by about 30% 
in the region of strange quark mass mga w 0.1, 
and y values of the Kyoto-Tsukuba and Kentucky 
groups become y = 0.46(14) and 0.42(3). This 
suggests a possible presence of large 0[mqa) ef- 
fects for the nucleon matrix element of strange 
quark density. 

3.3.3. baryon mass splitting and reduced 
matrix elements Fs and Ds 

For quenched QCD lattice estimates for the a 
term reasonably agree with the phenomenological 
value. Within the first-order flavor SU (3) break- 
ing formula , a large value of y is not consis- 
tent with this result unless the mass splitting of 
the baryon octet, or equivalently the matrix ele- 
ment {N\uu + dd-2'ss\N) « {N\uu + dd\N)conn., 
is small. 

To examine this point we recall the definition 
of reduced matrix elements Fg and Dg given by 



{B\q\^q\B) 



FsTr(i3t[A8,B]) 
DsTt{B^{\^,B}). 



(11) 



To first order in flavor S'L/(3) breaking, we have 
the relations Fs = (Ms — Mn) /2{ms — rh) and 
Ds = {Mb. + Mn ~2MT)/2{ms-m), from which 
we find experimentally that 
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Figure 7. Reduced scalar matrix elements Fs and 
Ds- For reference see Table 1. 

These values are to be compared with lattice 
estimates obtained from the relations, 



Fs 



{N\uu\N) 

conn. 

{N\uu-2dd\N)^ 



(13) 



F°"P- = 1.52, 



Df^- = -0.518. 



(12) 



where only connected contributions are kept since 
disconnected contributions are to be ignored to 
leading order of flavor symmetry breaking. 

Lattice results from the studies in Tab le p] are 
plotted in Fig. |^. Previous estimates pO|j21[ | are 
similar with results in this figure. While there 
may be some trend of increase toward the chiral 
limit, both matrix elements are small compared 
to experiments. It should be added that the ratio 
Fs/Ds is consistent, however. 

One can estimate the baryon mass splitting 
from results for the reduced matrix elements. At 
/3 = 5.7, for example, one finds Ms — = 
177(8)MeV and My. - Mn = 114(5)MeV |l|] 
as compared to the experimental values 379MeV 
and 254MeV. On the other hand, a recent spec- 
trum calculation aX (3 = 6.0 on a 32"^ x 64 
lattice reported Ms - Mn = 300(27)MeV and 
My-Mn = 185(17)MeV 111, which show a much 
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Figure 8. Connected contribution ASconn. = 
(Aii + Ad)conn. of the axial vector matrix element 
as a function of (m^/mp)^. For reference see Ta- 
ble 1. 

better agreement with experiments. 

Systematic analyses calculating both the scalar 
matrix elements and the baryon mass splitting 
within the same simulation have not been carried 
out so far. These are needed to elucidate the ori- 
gin of the problem in these quantities, especally 
to see whether it can be ascribed to a large vio- 
lation of scaling. 



4. NUCLEON AXIAL 
TRIX ELEMENTS 



VECTOR MA- 



4.1. Quark content of proton spin 

The forward matrix clement of the axial vector 
current 97^759 measures the contribution Aq of 
a given quark flavor q to the spin of the proton, 



{Ps\qi3i5q\Ps) =^ s- Aq, 



(14) 



where \Ps) is the proton state at rest with the spin 
projection in the z direction equal to s/2. Recent 
experiments on polarized deep inelastic lepton- 
nucleon scattering from SMC ^ and E143 [|8| 
have confirmed the earlier finding of EMC [59 
that the fraction of proton spin carried by quarks 
has a smaU value AE = Au-|-Ad-|-As « 0.2-0.3 
and that the strange quark contribution is non- 
vanishing and negative As w —0.1. These values 
are quite different from naive expectations from 
quark models. We tabulate the experimental val- 
ues and those from quark models in Table 3 (a) 
and (b). 
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Figure 9. Ratio R{t) for the (a) disconnected 
and (b) connected amplitudes of the axial vector 
matrix element Au evaluated by the wall source 
method without gauge fixing [ pTp^ for ii'=0.164 
at /3=5.7 on a 16"^ x 20 lattice. Solid lines are 
linear fits over 5 < t < 10. 

Recent lattice results for the connected contri- 
bution AEconn. — (Am -I- Ac?) conn, are compiled in 
Fig. ^ where we used the tadpole-improved renor- 
malization factor given by |32 



Za = [l- 



3K 



1 



(15) 



Lattice parameters for individual runs are sum- 
marized in Table 1. For comparison we also plot 
the experimental value of SMC at {m.,r/mpf' — 
0. The connected contribution alone is too large 
to explain the small value of AE. 

In Fig. H we show results for the R ratio 
for the up quark contribution from Rcf. |]ll| , |l3| 
where the disconnected contribution (Fig.^(a)) 
is obtained with the method of wall source with- 
out gauge fixing. Compared to the case of the 
scalar matrix element, the data for the discon- 
nected contribution are much noisier. The nega- 
tive value of the matrix element, not predictable 
in quark models, is clearly observed, however. A 
similar trend is also seen in the result obtained 
with the Z2 source shown in Fig. |l^ |^,|5|. 
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Figure 10. Ratio R{t) for the disconnected ampli- 
tude of the axial vector matrix element Au evalu- 
ated by the Z2 source method for if =0.154 
at (3=6.0 on a 16^ x 24 lattice. Solid line is a linear 
fit over 10 < i < 13. 

Table 3 

Axial vector matrix elements from (a) experi- 
ments, (b)quark models and (c) lattice calcula- 
tions. 







Au Ad 


As 


AS 


(a) Ej 

SMC 

E143 


cperiment 

33 0.80(6) -0.46(6) 
3|] 0.82(6) -0.44(6) 


-0.12(4) 
-0.10(4) 


0.22(10) 
0.27(10) 


(b) Quark model 
non-rel. 4/3 -1/3 
rel. 1.01 -0.25 






1 

0.76 


(c)La1 
Ref. 1 
Ref. 


ti 
13 
15 


ce results at = 

1 0.638(54) -0.347(46) -0.109(30) 0.18(10) 
1 0.79(11) -0.42(11) -0.12(1) 0.25(12) 



Values of the connected and disconnected con- 
tributions obtained in ref. |Q are plotted in 
Fig. ^ for u, d and s quarks separately as 
functions of {mT^/nip)'^. The magnitude of dis- 
connected contributions is small for each quark. 
However, the total amplitude AS receives a neg- 
ative disconnected contribution from u, d and 
s quarks. As a result the value of AS is sig- 
nificantly reduced from that obtained from con- 
nected contributions alone. This is shown in 
Fig. |lj where filled symbols are the total am- 
plitude and open ones connected contributions 
taken from Fig. |8[ We observe that the total am- 
plitude points toward the experimental value for 
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Figure 11. Axial- vector matrix elements for u, 
d and s quarks as functions of {mT^/mp)^ p^ . 
Strange quark mass is fixed to the physical value. 
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Figure 12. Total value of AS as a function of 
{m.^ /nipY . For reference see Table 1. Triangle is 
a result using the anomaly equation (5). 

light quarks. 

A more detailed comparison is made in Table 
3 where the total contribution extrapolated to 
niq = is listed for each quark. Lattice results 
are in reasonable agreement with those from ex- 
periments. 

4.2. Isovector axial coupling gA 

The SU{'i) reduced matrix elements for the ax- 
ial vector current are given by 



Fa = ^(A^)con„. 
Da = i(Au-2Ad)^ 



(16) 
(17) 
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Experimental values of these matrix elements are 
deduced using neutron /3 decay, which gives a very 
precise value of the isovector axial coupling gA, 



g''^^- = Fa + Da = 1.2573(28) 



(18) 



and hyperon (3 decay, which yields 
i^A - -Da -0.327(20). (19) 
From the two estimates one finds for the ratio, 
Fa/Da - 0.58(5) (20) 

In Fig. |l^ we plot recent lattice results for gA 
and the ratio Fa/ Da- We see that lattice data for 
gA are consistently smaller than the experimental 
value by about 20-30 % toward light quarks, while 
the ratio Fa/ Da is consistent. 

We remark that a previous result ^ for a \/3- 
blocked quark action shows a simiar trend if a 
tadpole-improved rcnormalization factor appro- 
priate for the action is employed. 

The disagreement of gA represents another 
problem with flavor non-singlet nucleon matrix 
elements. It does not seem to arise from the ef- 
fect of quenching since full QCD estimates plot- 
ted by upward triangles do not deviate from 
quenched results. Effects of scaling violations 
also seem small since results from two values of 
f3 at 5.7 and at 6.0 |l8|j3^j3l| are in agree- 
ment. One may suspect a possible ambiguity in 
the choice of the rcnormalization constant Za, 
especially because the ratio Fa/ Da agrees with 
experiments. However, the same rcnormalization 
factor enters in a calculation of the pion decay 
constant /jr, for which quenched lattice results 
a,t (3 — 5.7 — 6.0 agree with the experimental 
value within about 10%. Furthermore a non- 
perturbatively determined value of Za obtained 
a,t /3 — 6.0 and K— 0.1515 and 0.153 is quite con- 
sistent with the perturbative result ( p^ pof . 

5. HADRON SCATTERING LENGTH 

The theoretical basis for a lattice calculation of 
hadron scattering lengths is the formula proven 
by Liischer |^l|; it relates the s-wave scattering 
length uq to the energy level E of the lightest two- 
hadron state with vanishing spatial momentum in 
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Figure 13. Isovector axial coupling gA and 
Fa/ Da as functions of [raT^ /mpf . For the mean- 
ing of symbols see Fig. 8. 

a cubic box of length L according to 

„ 47ran , an fln , „ , 1 x / x 

i.-2™.--^(l + c,^ + c.^)-fO(^)(21) 

with ci = -2.837297 and C2 = 6.375183. The 
technical difficulty of calculating general four- 
point functions (including disconnected ones) 
needed for an extraction of the energy level E 
was overcome by the wall source method without 
gauge fixing, and tt-tt, tt-TV, K-N and N-N scat- 
tering lengths have been calculated in quenched 
QCD While no new simulations have been 

carried out since LATTICE 93, two theoretical 
studies aiming at a deeper understanding of the 
lattice calculation of hadron scattering lengths 
have been made (see also Refs. [^,|3| for related 
investigations). 

5.1. TT-TT scattering length 

Chiral perturbation theory suggests that quan- 
tities calculated in quenched QCD often develop 
spurious infrared divergences in the chiral limit, 
which are not present in full QCD jij]. Bernard 
and Golterman investigated the effect of 
quenching in the energy level E of the two-pion 
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Figure 14. N-N scattering length oq calculated 
by a one-boson exchange model for ^Sq and ^Si 
channels as functions of quark mass niq. The 
physical point corresponds to = 4.8 MeV. 



state using quenched chiral perturbation theory 
technique, extending a previous work on the sub- 
ject They found that enhanced finite- volume 
corrections of order L° = 1 and L^^ occur in E 
at one loop due to the double pole in the rj' prop- 
agator. This implies that it is not consistent to 
use Liischer's formula in quenched QCD. 

On the other hand, a quenched lattice calcu- 
lation of TT-TT scattering lengths with the Kogut- 
Susskind quark action |^ yielded results consis- 
tent with the prediction of current algebra and 
PCAC, which is known to hold for the Kogut- 
Susskind case |46). They suggest that this agree- 
ment is due to an accidental choice of lattice pa- 
rameters such that the enhanced finite-volume 
corrections are not large for the values of m^L 
of the simulation. 

5.2. N-N scattering length 

The nucleon-nucleon scattering length is a 
fundamental parameter characterizing the low- 
energy properties of the nuclear force. Experi- 
mentally they are well determined and are known 
to have large values. 



ao(^'5'o) 



-20.1(4) fm, 
-5.432(5) fm. 



(22) 



esting dynamical problem. However, a number of 
obstacles exist toward realistic calculations. A 
very large lattice is needed to apply Liischer's 
formula in order to suppress 0{L~^) corrections 
neglected in ([2l|). Extraction of the scattering 
length in the '^Si channel requires a calculation 
of the lowest scattering state orthogonal to the 
ground state, which is the deuteron bound state. 

A possible strategy in this situation would be to 
examine the quark mass dependence of the scat- 
tering lengths starting from a heavy quark where 
simulations are much easier. For this approach 
we need to understand the behavior of scattering 
lengths for heavy quarks. 

For heavy quarks the N-N interaction becomes 
shorter ranged since pious exchanged between nu- 
cleons are heavier. This will force the nucleons 
out of the attractive well of the potential, which 
suggests that the deuteron becomes unbound as 
the quark mass increases. 

A closer examination of this possibility was 
made in Ref. taking a phenomenological 

model of one-boson exchange potentials and 
varying hadron masses in the model as a func- 
tion of quark mass according to lattice data. The 
results for the N-N scattering lengths are plot- 
ted in Fig. |lj. The divergence of the scattering 
length for the '^Si channel taking place at niq = 
6.3 MeV signals unbinding of the deuteron. For 
TTiq > 6.3 MeV, both 00(^6*0) and ao('^S'i) are pos- 
itive and rapidly decrease to have values of order 
1 fm. Also the value of ao(^S'i) is larger than that 
for aQ{^So) for heavy quark indicating a stronger 
attraction in the ^Si channel. 

These are precisely the features found in a lat- 
tice calculation of N-N scattering lengths 
The results are shown in Fig. 15 together with 



Since there are no constraints from chiral symme- 
try for N-N scattering, explaining these values 
from first principles of lattice QCD is an inter- 



those for n-n and tt-N cases for comparison. 

6. CONCLUSION 

In the last two years, a substantial progress 
has been made in the lattice QCD study of the 
flavor singlet nucleon matrix elements. Appli- 
cations of efficient calculational techniques for 
handling disconnected quark loop amplitudes re- 
vealed significant sea quark contributions in the 
flavor singlet scalar and axial vector matrix el- 
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Figure 15. N-N scattering length in physical 
units at (3=5.7 on a 20^ lattice. Also shown are 
the TT-TT and n-N scattering lengths at the same 
/3 on an 12^ x 20 lattice. 

ements. Adding sea and valence quark contri- 
butions the total magnitude of these matrix ele- 
ments are comparable to experimental estimates. 
Although the nature of systematic errors such as 
scaling violations and quenching effects need to 
be understood more precisely, the encouraging re- 
sults found so far point toward an eventual reso- 
lution of the problems related to the flavor singlet 
matrix elements within lattice QCD. 

In the course of the studies it has become ap- 
parent that current lattice results for flavor non- 
singlct scalar and axial vector matrix elements 
differ from experiments beyond statistical errors. 
A systematic scaling study appears needed to pin- 
point the origin of the problem. 

We have also discussed some recent work in the 
study of hadron scattering lengths. The model 
result that the deuteron becomes unbound for a 
quark mass only slightly larger than the physical 
value suggests that much interesting physics can 
be done on the subject of N-N scattering. At the 
same time the possibility of a severe quenching 
error pointed out for the tt-tt case signifies the 
need of caution in the use of quenched QCD. 

Study of flavor singlet physics is a young sub- 
ject in lattice QCD. Results obtained so far 
should serve as a stepping stone for further de- 
velopment in this area in the years to come. 
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